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QCD factorization with heavy quarks 
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We further analyze the definition and the calculation of the heavy quark impact factor at next-to-leading (NL) 
log s level, and we provide its analytical expression in a previously proposed fc-factorization scheme. Our results 
indicate that fc-factorization holds at NL level with a properly chosen energy scale, and with the same gluonic 
Green's function previously found in the massless probe case. 



1. INTRODUCTION 

Recent improvements Q| of the next-to-leading 
logs (NLx) results Q in the BFKL framework, 
have stabilized the small-x behaviour in QCD, so 
that a phenomenological analysis of deep inelastic 
processes (DIS) seems now possible. 

However, both the gluon density (satisfying 
the improved equation) and the impact factors 
are needed in order to use fc-factorization || 
to compute DIS or double DIS processes. So 
far, NLa; impact factors have been found for the 
unphysical case of massless initial quarks and 
gluons only f^,||. Partial features for massive 
quarks |||| and for colourless sources (?J are 
known too. 

In this talk we present the complete results 
for the case of initial massive quarks derived in 
Ref. § . These results allowed us to check the va- 
lidity of the fc-factorization scheme introduced in 
Ref. |J], or, in other words, to derive probe inde- 
pendent gluon Green's function with an explicit 
massive quark impact factor which satisfies the 
expected collinear properties. Furthermore, we 
developed as a byproduct some analytical tech- 
niques which are needed to deal with two-scale 
problems, which are hopefully useful to cope with 
the physical cases also. 
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2. k- FACTORIZATION IN DUET PRO- 
DUCTION 

Following ||, the colour averaged differential 
cross section for the high-energy scattering of two 
partons a and b is factorized in a gauge-invariant 
way into a Green's function Q u and impact factors 
h a and h b (Fig. 0) 



d fl- 



ab 



d[fci]d[fc 2 



duj 

2niui \so(ki,k 2 ) 
x h a {k 1 )g u (k 1 ,k 2 )h b (k 2 ) 



(1) 



The transverse momenta fci and k 2l defined with 
respect to the incoming momenta p\ and p 2l play 
the role of hard scales of the process. 

At the next-to-leading log a; (NLx) accuracy 
the Green's function Q u has the following general 
form 



Q u = (l + a s H L ) 



x (1 + a s H R ) , 



(2) 



where K$ and Kn l are the leading log x (La;) and 
the NLz BFKL kernels §] respectively, H R (H L ) 
are operator factors introduced in [Q so as to pro- 
vide partonic impact factors free of double log 
collinear divergences and a s = a s N c /Tr is the di- 
mensionless strong coupling constant. 

As explained in j| , the identification of the sec- 
ond order impact factors, ha and h b , is affected 
by a double factorization scheme ambiguity, due 
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Figure 1. Diagrammatic representation of k- 
factorization. 




Figure 2. Real gluon emission in the fragmenta- 
tion region of quark a . 



to both the choice of the scale sq and of the ker- 
nels H R (H L ). 

3. FACTORIZATION SCHEME AND 
CALCULATIONAL PROCEDURE 

Let's consider first the high-energy scattering 
of two partons a and b where a = q is a heavy 
quark of mass to with real emission of an extra 
gluon g that we assume in the heavy quark frag- 
mentation region (Fig. |2|). In terms of invariants 
S2 = (q + Pa) 2 > si = (p 3 + q f. The Born dif- 
ferential cross section in this high energy region 



was calculated in [g . Though complicated at first 
sight it reduces, as expected to the known ||J re- 
sult for to — ► 0, and matches the Lx differential 
cross section 



da 



qgb 



dz\ d[fei] d[fc 2 



= / l (°)( fel )/ l (0)( fe2 ) 
1 



«r»r(i 



s)fl 2e Z\ 



- , (3) 



in the limit z\ — > 0, being z\ the momentum frac- 
tion of k\ = pi — P3 with respect to the incoming 
momentum p±, M°)(fe,) the leading order impact 
factor and q = fci + 

However, as pointed out in ||, for eq.(||) to be a 
good approximation to the total result, we should 
require (q = \q\,h = \k t \) 



z\ <C q/ki , q/m , ki/m 



(4) 



The first two cutoffs can be summarized by Z\ < 
q/ Alax{k\ 1 to), which is a coherence condition for 
the case of heavy quarks, saying that the rapidity 
of the gluon cannot exceed that of the final quark. 

By integrating the leading expression (^) with 
the constraints (^) in the fragmentation region 
z\ > q/\fs, an estimate of the leading contri- 
bution contained in the complete result, which 
should be subtracted out in order to yield the 
impact factor in the massive quark case, was 
found ||. Then, by considering both real and 
virtual contributions to the fragmentation func- 
tion Fq(z\,h\,k2), we introduced the following 
definition of the impact factor hq (k): 



/ dzi / d[fci]F q (zi,fci,fe 2 ) 



h^(k 2 ) + / d[k 1 ]a s hW{k 1 )K {k 1 ,k 2 ) 



log- 



log — e qkl 

ki 



(5) 



Max{k\, to) 

Compared to the subtraction (or factoriza- 
tion) scheme adopted in Q] for m = 0, the ex- 
pression (H) differs by the replacement k\ —> 
Max(ki, to), which leads, by adding the symmet- 
rical fragmentation region, to the choice for the 
factorized scale in cq. (|lj) 



sq = Max{ki, mi)Max(k2, TO2) 



(G) 
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mi being the mass of quark a and m 2 the mass of 
quark b, and in particular contains the subtrac- 
tion term log q/kiQ q k t which provides the expres- 
sion (H R = hI=h) 

for the H kernel in the fc-factorization formula. 

In order to simplify the subsequent calcula- 
tions, the known result for m — was used 
and only the difference for a non vanishing mass 

AF q (z 1 ,fc 1 ,fe 2 ) = 

= F q (* 1 ,fc 1 ,fc 2 )-F q m=0 (* 1 ,fc 1 ,fc 2 ) , (8) 

was explicitly computed. Then, we found the fol- 
lowing relationship between the massless quark 
and the heavy quark impact factors 

h£Hk 2 ) = h { H =0 (k 2 ) + (9) 
+ J dzij d[k 1 ]AF ct (zi,k 1 ,k 2 ) 

+ J d[fci]5, ^ q 0) (fci)^o(fci,fe 2 ) log^9 mfcl 

Notice that the integration limits in z\ have been 
extended down to z\ — 0. Since AF q is regular 
at Z\ = this change introduces only a negligible 
error of order 1/s. 

4. MELLIN TRANSFORM AND ITS IN- 
VERSE 

In order to perform the calculation outlined 
in eq.(^), we proceeded in two steps. First, the fei 
integration was performed analytically by reduc- 
ing the fei-integrals to two denominators. Then, 
the virtual contribution || was considered and or- 
ganized in terms of momentum fraction integrals 
only. Finally summing up real and virtual contri- 
butions to the fragmentation vertex we obtained 
an expression for the difference AF q (fc2), aris- 
ing from the second term in the r.h.s. of eq.(^|). 
To perform the last integrations we calculated its 
Mcllin transform 

AF q ( 7 ) =T(l+e) (m 2 )~ £ 



/ d[fc2] (5) 7 AFq(fc 



which allowed us to disentangle the (m/ In- 
dependence, yielding 

AF q ( 7 ) = A £ (m 2 Y 

r( 7 + e)r(l - 7 - 2e)r 2 (l - 7 - e) 
X 8r(2 - 2 7 - 2s) 

1+e 2 



7 + 2e 1 - 27 - 4e 

1 1 



1 - 7 - 2e 3 - 2 7 - 2e 



, (10) 



where A £ is a constant that contains the depen- 
dence on the strong coupling constant and some 
colour factors. 

It is straightforward, though not trivial, to 
show that eq. ([!(]) converges only in the small 
band 1 — 2e < Rej < 1 — e. The inverse Mellin 
transform was thus defined as 



AF q (fc 2 ) = \\ 
m z J 1 



dj 
2^i 



-2e<fle7<l-e 
2 \ -7-e 
\ I — ?r I AF p ( 7 ) . 



in 



Then, displacing the integration contour around 
the positive or the negative real semiaxis, i.e. en- 
closing all the poles placed either at 7 > 1 — e 
or 7 < 1 — 2e, we calculated the different cor- 
rections of order 0(m/k 2 ) n or 0{k 2 /m) n to the 
impact factor in the limits k 2 > m 2 and k 2 < m 2 
respectively. 



5. IMPACT FACTOR 

Our final result for the heavy quark impact fac- 
tor at the next-to-leading level reads 

h q (k 2 )^h^(k 2 )\ smg + h q (k 2 )\ fmtte , (11) 
where the singular piece is defined as 



(12) 



h^{k 2 



a s N c 
2tt 



4 l0 S~! ) Qk 2 m 

2 m z 
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and 



h q (k 2 



K 



I finite 9 

tt2 / :-! 
6 



2tt 



Res[h( 7 )})e k 



Re~f>l 



2+ Y, Res[Mi)])e 

Rey<l 



m k -2 



(13) 



is the finite contribution, with 



K = --- 

18 6 9N, 



%L . (14) 



As for the massless case, the singularities propor- 
tional to (11/6 — n//3iV c ), the beta function, were 
absorbed by the running strong coupling constant 
a s (/c2). The function ^1(7) provides the correc- 
tions II of order 0(m/fc 2 ) and 0(fc 2 /r7i) to the 
impact factor for k\ > m? and k\ < m 2 respec- 
tively. 

Notice that all double log contributions of type 
1/e 2 and 1/ e\og{k\/m 2 ) appearing in the inter- 
mediate steeps of the calculation canceled out 
which means that indeed our subtraction of the 
leading kernel was effective, thus lending credit 
to the scale (|J) and to the H- kernel (0). The 
remaining singularities of the impact factor are 
single logarithmic ones ~ 1/e. In fact, the 
impact factor is actually finite, with the ex- 
pected \og(k 2 /m 2 ) dependence predicted by the 
DGLAP equations, the divergent piece 5/4 (£2) 
in eq . (|l2l) , see M, can be interpreted as a finite 
mass scale change, i.e. the scale leading to a finite 
massive quark impact factor differs from eq. @ by 
a finite renormalization of the quark mass, which 
is a normal ambiguity in this type of problems. 

6. CONCLUSIONS 

Starting from the explicit squared matrix ele- 
ment for gluon emission we motivated the sub- 
traction of the leading term, and we performed 
the k\ and z\ integrals needed to provide an ex- 
plicit result for the heavy quark impact factor. 

Even if the cross section being investigated is 
unphysical, the relevance of our results stems 
from the consistency of the following features: 



(i) the validity of the fc-factorization formula (|l|) 
with scale so = Max(ki, mi)Max(k2, TO2); (ii) 
the explicit expression of the impact factor with 
factorizable single logarithmic collinear diver- 
gences, and (iii) the probe-independence of the 
subleading H- kernels of the CC scheme Q, de- 
fined in eq.(^|). 

Of course, the real problem is to provide an ex- 
plicit expression for the DIS impact factors. But 
- if the lesson learned form the L and NL ker- 
nels is still valid - the impact factor's magnitude 
is not expected to be much different from their 
approximate collinear evaluation. 
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